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^ Abstract 

It is shown that the data on the pion charge form factor admit the possibihty for a substantial 
sea-quark components in the pion wave function. If the charge form factor is calculated with 
instant form kinematics in a constituent quark model that is extended to include explicit {qqf' 

^ : 

components in the pion wave function, that component will give the dominant contribution to the 
calculated vr^ charge form factor at large values of momentum transfer. The present experimental 
values can be described well with {qq)"^ component admixtures of up to 50%. The sensitivity 
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' of the calculated vr"^ charge form factor to whether one of the quarks or one of the antiquarks is 
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taken to be in the P-state is small. 
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I. INTRODUCTION 



The small mass of the pion, in comparison to that of the vector mesons, suggests that 
the pion is a collective state of (qq)^ configurations, with many values of n. Therefore it 
is something of a riddle that it is possible to describe the empirical charge form factors of 
the charged pions satisfactorily with phenornenological wave functions under the assumption 
that they are pure quark-antiquark states [1]. It is then a natural question to ask whether 
the empirical pion charge form factors can exclude the presence of the expected multiquark 
configurations in the pion. 

The negative parity of the pion requires that in the simplest sea-quark configuration, 
[qqY, at most 3 of the constituents can be in the ground state and that either one of the 
quarks or one of the antiquarks is raised to the P— state (or a higher odd— L— state). As 
this is energetically unfavorable it suggests that the probability of that configuration may 
be small in comparison to that of the qq component. The situation is analogous to that 
in baryons, for which positive parity requires that in a qqqqq admixture either one of the 
quarks or the antiquark has to be in the P— state I2I. 

Here the charge form factor of the (charged) pion is calculated in an extension of the 
constituent quark model to include admixtures of the simplest sea-quark configurations 
(gg)^ in instant form kinematics. The calculation is made both for the case where one of 
the quarks and where one of the antiquarks is in the P— state. It is in both cases found 
that inclusion of the sea-quark configuration allows a good description of the empirical form 
factor, even if it represents as much as half of the wave function. The main point is however, 
that as soon as there is a non vanishing probability for the sea quark component in the wave 
function, that component will lead to the dominant contribution to the charge form factor at 
large values of momentum transfer in the case of instant form kinematics. As a consequence 
the fact that it is possible to achieve a quantitatively satisfactory fit to the empirical charge 
form factor with wave function model for the conventional qq component alone does not rule 
out the presence of significant sea-quark components in the pion. 

The method of calculating the charge form factor in the constituent quark model with 
instant form kinematics developed here can readily be extended to sea-quark configurations 
with larger numbers of sea-quark qq components. 

The configurations of the [qqY system that are possible in the pion are described in 
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TABLE I: The flavor-spin-color configurations of the qq and qq pairs. 













qq 


SU ii) flavor 


[2]f, [11]f 


[2]f, 


SU (i)gpin 


[2]5, [n]s 


[2]5, [n]s 


SU {2i)color 


[2]c, [ll]c 


[22]c, [211]c 



section 2. Section 3 contains a description of the pion wave function and the different form 
factor contributions. The calculated results for the pion charge form factor are given in 
section 4. Finally section 5 contains a summarizing discussion. 

II. LIGHT FLAVOR {qqf CONFIGURATIONS IN THE PION 

The lightest {qgY component in the pion meson contains only the light flavor quarks 
u and d, which form the fundamental representation of the SU{2) flavor symmetry. The 
flavor-spin-color conflgurations of the qq and qq are listed in tabled by their Young patterns. 

The wave functions of the qq and qq pairs in the flavor-spin-color-orbital space should 
be totally antisymmetrized, respectively. In addition, the odd parity of the pion meson 
requires that either a quark or a antiquark in the qqqq component is in P-state (or higher 
odd— L— state. 

In the case where one quark is in the P-state and both antiquarks are in the S-state the 
flavor-spin-color wave functions of the qq pairs in the {qq)'^ component are totally symmetric 
while that of the qq pairs are totally antisymmetric. This leads to four possible color singlet 
conflgurations of the (qq)"^ component with = 0~: 

a. {[2]^[2]5[2]c},,{[2]i.[ll]5[22]c},, 

b. {[2]p[2]s[2]cU{[U]f[2]s[22]c},, 

c. {[2UU]s[n]c}q,mF[2]s[2n]ch, 

d. {[n]F[2]s[n]c}MF[2U2n]c},,. (i) 

It is natural to assume that the (qq)"^ conflguration with the lowest energy shall have the 
largest probability in the pion besides that of the conventional qq component. The splitting 
of the energy of the (qq)"^ components ([T]) is determined by the hyperflne interaction between 
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TABLE II: The expectation values —jr{a\Hj\a) of the {qq)"^ configurations in eq. ([T]). 



a 


a 


b 


c 


d 


— ^{a\Hj\a) 


8 
3 


16 


16 
3 






the quarks and the antiquarks. This interaction will here be taken to have the schematic 
form: 

HJ = -XY.^^■^A'i■^■ (2) 

Here a are the Pauli matrices, A*^ are the Gell-Mann matrices in color space and X is a 
positive constant with the dimension of energy. This interaction has the same color and 
spin dependence as the color magnetic hyperfine interaction, which arises from single gluon 
exchange. The contributions to the energies of the 4 configurations ([T]) that arise from the 
schematic hyperfine interaction ([2]) can be determined by the re-coupling method described 
in [sl and 4| and are listed in table HTl 

In the case where one antiquark is in the P— state, and both the quarks are in the iS— state 
the wave functions of the qq pairs in fiavor-spin-color space have to be totally antisymmetric, 
while those of the qq pairs shall be totally symmetric. The possible {qq)^ configurations in 
the pion are in this case the following: 

a'. {[2]F[2]s[n]c}MF[n]s[2n]ch, 

b'. {[2]^[2]5[11]c}..{[11]f[2]5[211]c},-, 

c'. {[2Ull]s[2]cUmF[2]s[22]ch, 

d'. {[U]p[2]s[2]cU{[2]F[2]s[22]ch, . (3) 

The expectation values of the hyperfine interaction between quarks in the {qq^ configura- 
tions in ([3]) are listed in table IIIIl The results in the tables HI] and IIIII show that the hyperfine 
interaction between quarks leads to the same energy levels of the [qqY configurations, in 
which the antiquark is in the S-state and P-state. This is a consequence of the fact that the 
hyperfine interaction is independent of the angular momentum of the constituent quarks. In 
the case, where the antiquark is in the P-state the [qq^ configuration d' has the lowest en- 
ergy, which is equal to that of the lowest energy {qqy configuration b, in which the antiquark 
is in the S-state. These two configurations are thus likely to constitute the most probable 



TABLE III: The expectation values —-^{a'\Hi\a') of the {qqY configurations in eq. ([3]). 
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-^{a'\Hi\a') 


16 
3 





8 
■3 


16 



{qqY components and therefore to be most significant in the structure of the pions. The 
roles of such configurations in the pion charge form factor will be considered in the following 
section. 

III. 7r+ CHARGE FORM FACTOR 

A. The form factor contribution from the qq component 

The present empirical results for the charge form factor of the vr"*" may be well described 
as a qq system in a Poincare covariant constituent quark model with the following orbital 
wave function model py]: 

0(^1, k2) = Af2, . (4) 

(1 + %^)'^ 

Here A/2g is a normalization factor, a and b are parameters and ki, i=l, 2, are the quark 
momenta in the rest frame of vr meson ki=0) . This wave function model will be 

adopted here for the qq component of the pion wave function. 

In the impulse approximation the contribution of the qq component to the pion charge 
form factor is obtained as the matrix element of the electric current density operator between 
the initial and final states in the Breit frames as: 

F(^^")(Q2) = J dV2\/«^e(Pl,Pl')0(^l,^2)</>(^l,fc2)- (5) 

The initial and final momenta of the constituents in their respective Breit frames are denoted 
Pi and Pi ' respectively {pi' = Pi + Q and p2 = ')• 

If the momentum transfer Q is taken to define the axis, the relations in instant form 
kinematics between the momenta in the Breit frames and the rest frames are: 

Pi± = ki_L = ki±' = pi± ' , 
Pi\\ = vokiii + v\i uji , 
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Pi\\ = ^okw + uJi , 
Ei = v\\ki\\ + VoUi , 

E'i = v\\k'i\\ + VQu'i . (6) 



Here the energy components are defined as 



uji = \l ki"^ + m? , = y ki"^ + m'^ , 



Ei = ^Pi'^ + m?, = \J pi"^ + m? . (7) 

In these relations m denotes the constituent mass and v = {vq, 6±, v\\} and v' = {v'q, 6±, Vy} 
the constituent boost velocities in the initial and final states. These satisfy the constraint 
= v"^ = —1. 

In instant form kinematics the boost velocities may be defined as : 

Q 



Q 



(8) 



Here n represents the number of constituents. 

The Jacobian that is induced by the transformation between the rest frame and the Breit 
frame of the meson is in the case of the qq component obtained as Q]: 

* = f(l-ll^). (9) 

The expression for the corresponding final state Jacobian J2 ©is obtained by replacement 
of the arguments by the corresponding primed coordinates. 
The electric current density operator in eq. is 



Se{P,P, 



l [E' + m){E + m) f p' ■ p 

M ^ AM^\ AE'E I {E' + m){E + m)^ ' ^ ' 



The expressions for the boost velocities of the constituents ([8]) reveal, that their magni- 
tudes fall with increasing number of constituents n, if the constituent mass is constant. Given 
that form factors fall with increasing momentum transfer it follows that, at sufficiently 
large Q^, the wave function component that contains the largest number of constituents 
will give the largest contribution to the form factor. This feature is explicit in instant form 
kinematics. It has a natural physical interpretation in that the form factor describes the 



probability that the system stays bound upon absorption of the momentum transfer Q. The 
relative probability for this to happen is smaller if few constituents absorb the momentum 
transfer than if many can share it so that the fractional momentum transfer per constituent 
is smaller. 



B. The form factor contribution from the (gg)^ component 

In the case where both of the antiquarks in the {qqY component are in the S"— state, the 
wave function of the {qqf' component, which has the symmetry configuration h in eq. ([1]), 
and which is expected to have the lowest energy, may be expressed as: 

= uu{du - ud) {[Is ® l's\i ® lx}o- {6c ® 6c}ic h, h) ■ (11) 

Here ki, i=l..A, are the momenta of the constituent quarks in the rest frame of the pion 
(J2i=i ki = 0). The spin triplet combinations of the qq and qq pairs are denoted I5 and I5, 
respectively. These combine with the P-state qq pairs to the total quantum numbers of pion 
= 0~ . The spherical harmonic for the qq pairs in the (qq)"^ component is defined as 

lx,m=^i„,, 6 = ;^(A;i - A-a) , (12) 

where (m = —1,0, 1) are the spherical components of ^i. In (fTTj) we have denoted the 
Young pattern representations of the color states of qq and qq pairs in table [T] with their 
corresponding dimensions (6). 

The orbital wave function of the [qq^ component is taken to have the form: 

$(^1,^2,4,^4) =Ar4, , (13) 

(1 + %^)^+^ 

where Af^q is a normalization factor. 

Then the contribution to the vr"'" charge form factor from the qqqq component in the case 
where one of the antiquarks is in the S-state then may be written as: 

F^MQ') = U'{Q')-^-B'iQ'). (14) 

Here the terms A and B are defined as: 

A'm = \j d%d%d%^Mi)Ji{i)Se{pi,if,)(i ■ ^Mh, h, h, hmk[, k'„ k'„ k'^, (15) 

B'iQ') = \j d%d%d%^M4)J!,{4)Se{pAJdii ■ e>(^i, h, h, hmk[, k'^, k'„ k'^).{16) 



Here A^{Q'^) represents the matrix element where a photon couples to the 1-st quark in the 
{qqY component, while B^[Q'^) represents the matrix element where the photon couples the 
antiquark (the fourth constituent). 

The Jacobians for the transformations between the corresponding Breit frames and the 
rest frames are: 

J.(4) = -«(l-„,^). (18) 

As in the case of eq. ([5]), the primed variables in eqs. (ITSl) and (fT6l) represent the final 
states variables that correspond to the initial state variables without primes. Here we do 
not take into account the Wigner rotation of the spin axis that is caused by the boosts, 
as its consequences are numerically insignificant for momentum transfers below 10 GeV^ 
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Comparison of the expressions for the Jacobians ([9]) and (jlSll for the case of 2 and 4 
constituents, respectively, makes it clear how to generalize these expressions to the case of 
n constituents. 

In the case where one antiquark is in the P— state and both quarks are in the ground state, 
the wave function of the {qqY component with the lowest energy, which has the symmetry 
configuration d! in eq. (I3l) may be obtained from that in the case where one of the quarks 
are in the P— state ( ITTl) by the replacements: 

u —d d u . (19) 

The explicit expression is then 

k+)p = -^{ud - du)dd{[ls ® l^i ® I'Jo- {6c ® 6c}ic Hh,k2, h, h) . (20) 
Here the spherical harmonic for the qq pair is denoted as 

l'x,m=e3„, ^3 = ^(4-4). (21) 

The explicit expression for the ir'^ charge form factor in the case where one antiquark is 
in the P— state are is: 

F^'4Q') = l^^iQ') + Is'm ■ (22) 
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Here the orbital integrals are defined as: 

A'^m = \j d%d%d%^Mi)J[ii)SM,Pi)C3 ■ CsHh, h, h, hMK, 4, K), (23) 

The symmetrical form of the expressions (fTT!) and (!20l) has the consequence that 

A^(g2) = B^{Q^) , BP{Q^) = A^{Q^) . (25) 

Because of the symmetry structure of the spin-flavor-color state of the {qqf' components, 
there is no contribution of off-diagonal {qqY Qq transition matrix elements to the pion 
charge form factor. In the case of the nucleons the contribution of such transition matrix 
elements to the form factors are much larger than that of the corresponding diagonal matrix 



IV. RESULTS 

To investigate possible role of the [qqY component in the form factor of the tt^ meson, 
the wave function parameters are chosen so that the combined contribution of the qq and 
the [qq)"^ components yield a form factor that agrees with the empirical one under the 
assumption of a probability for the latter component of 10 %. The corresponding wave 
function parameters a, 6 (jlj) and A, B f|T3l) are listed in table [IVl The constituent quark 
mass was taken to be 120 MeV, which is close to the value required to describe the nucleon 



The 



form factors in instant form kinematics with a wave function of corresponding form 
wave function parameters for the qq and the {qqf' components were chosen such, that the 
empirical mean square radius for the pion, = 0.44 fm^, was recovered. 

The calculated result for the vr"*" charge form factor is shown in Fig. [1] for the case where 
both antiquarks are in the iS— state in the {qqf component. The result indicates that above 
1 GeV^ that with these parameters the main form factor contribution arises from the smaller 
{qcif component. That this should be so is in fact quite natural, as in the case of elastic 
form factors, the form factor falloff with momentum should depend on divided by the 
square of the number of involved constituents. In this case the contribution of the {qqY 
component is very small (and in fact negative): -0.03 fm^. The sign of this contribution 
depends on the parameter values. 
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FIG. 1: The 7r+ charge form factor obtained with a 10% (qq)'^ component probabihty with both 
antiquarks in the S— state. Sohd hne: The contribution from the qq component; Dashed hne: The 
contribution from the {qq)'^ component; Dotted hne: the result combining the contributions from 
the qq and the {qq)"^ components. The data sets CERN 1986, DESY 1978, JLAB 2006 and JLAB 



2007 are taken from refs. 



Id, 



Hi, respectively. The data point at = 1.60GeV^ from ref. 



Ill ] has been shifted for better visibility. 



The corresponding calculated results for the vr^ charge form factor for the case, where 
one of the antiquarks is in the P— state are shown in Fig. O These results are in fact very 
similar to those obtained in the former case, where both antiquarks are in the S*— state, the 
main difference being a slightly larger magnitude for the contribution to the mean square 
radius from the (g^)^ component in the present case (-0.04 fm^). 

In Fig. [3]the form factor is shown as obtained for different values of the probability for the 
(qq)'^ component. These results were obtained by only slight variation of the 2 parameters 
in the wave function of the (qq)"^ component (|T3l) . These results show that the present 
empirical data on the pion charge form factor can allow for a {qqY component probability 
of up to 50%. 
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TABLE IV: The model parameters. 



mq (MeV) 



b (MeV) 



B (MeV) 



A 



10% 
20% 
40% 
50% 



120 
120 
120 
120 



190 
190 
190 
190 



100 
110 
139 
143 
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2.3 
2.3 
2.3 
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FIG. 2: The vr^ charge form factor with a 10% probabihty for the {qq)'^ component with one 
antiquark in the P-state. The labehng of the curves is the same as that in Fig. [H 



For comparison the form factor that is obtained for the pure qq quark model for the pion 
is shown in Fig. HI These resuhs were obtained with the constituent quark mass value 80 
MeV and with the parameters a = 2.0 and b = 198 MeV in the qq wave function model (jlj). 
This shows that the pion charge form factor may be described with such a simple model 
wave function in instant form kinematics, a result that was noted for the case of front form 
kinematics in ref.jl|. 
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FIG. 3: The vr^ charge form factor obtained with a 10%-50% (qq)'^ component probabiHties with 



both antiquarks in the S'-state. The data sets CERN 1986, DESY 1978, JLAB 2006 and JLAB 

y, 3, 



2007 are taken from refs. 



Ill ], respectively. The data point at Q = 1.60GeV from ref. 



Ill ] has been shifted for better visibihty. 



V. CONCLUSIONS 

The results of this study shows that the present data for the charge form factor of the 
charged pions may described as well with inclusion of sea-quark configurations with probabil- 
ities up to least 50 % in the covariant constituent quark model with instant form kinematics 
as without such configurations. 

When the sea-quark component is included in the form factor, that component will give 
the dominant contribution to the form factor at sufficiently large Q^, which dominates over 
that from the qq component, however small its probability. The [qq^ component corresponds 
to structures that have shorter range than the basic qq component. This is illustrated in 
FiglH where the charge density contributions from the qq and the {qq^ components, along 
with their sum, is show for the case, in which the probability of the latter component is 
20%. If the probability of the (qq)^ component is taken to be larger, the peak in the profile 
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FIG. 4: The it~^ charge form factor for the pure qq model. The labehng of the curves is the same 
as that in Fig. [H 



r'^p{r) moves towards that of the qq component. 

The fact that the sea-quark contribution gives the largest contribution to the form factor 
at large values of momentum transfer is a consequence of the fact that the magnitude of the 
energy denominator in the expressions for the boost velocities for the constituents ([8]) grows 
with the number of constituents. This implies that the momentum transfer is shared by the 
largest number of constituents, and in effect, to a smaller relative momentum transfer per 
constituent. As the form factor is a monotonically falling function of momentum transfer, 
the consequence is that the largest contribution at large is given by the component with 
the largest number of constituents. 
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FIG. 5: The "charge density" of the 7r+ for a (qq)^ component with 20% probabihty. The contri- 
butions from the qq and {qq)^ components are denoted 2Q and AQ respectively. 

[4] R. L. Jaffe, Phys. Rev. D15, 267 (1977); R. L. JafFe, Phys. Rev. D15, 281 (1977) 

[5] B. Juha-Di'az, D.O. Riska, F. Coester, Phys.Rev. C69, 035212 (2004) 

[6] F. Coester and D. O. Riska, Nucl. Phys. A728, 439 (2003) 

[7] Q. B. Li, and D. O. Riska, Nucl. Phys. A791, 406 (2007) 

[8] S. R. Amendolia et al., Nucl. Phys. B277, 168 (1986) 

[9] H. Ackermann et al., Nucl. Phys. B137, 294 (1978) 

[10] T. Horn et al., Phys. Rev. Lett. 97, 192001 (2006) 

[11] V. Tadevosyan et al., Phys. Rev. C75, 055205 (2007) 



14 



